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Abstract 



We calculate the critical value of the hopping parameter, n c , in Lattice 
QCD, up to two loops in perturbation theory. We employ the Sheikholeslami- 
Wohlert (clover) improved action for fermions and the Symanzik improved 
gluon action with 4- and 6- link loops. 

The quantity which we study is a typical case of a vacuum expectation 
value resulting in an additive renormalization; as such, it is characterized by 
a power (linear) divergence in the lattice spacing, and its calculation lies at 
the limits of applicability of perturbation theory. 

Our results are polynomial in csw (clover parameter) and cover a wide 
range of values for the Symanzik coefficients q. The dependence on the num- 
ber of colors N and the number of fermion flavors Nf is shown explicitly. In 
order to compare our results to non perturbative evaluations of k c coming 
from Monte Carlo simulations, we employ an improved perturbation theory 
method for improved actions. 

Keywords: Lattice QCD, Lattice perturbation theory, Hopping param- 
eter, Improved actions. 

PACS numbers: 11.15.Ha, 12.38. Gc, ll.10.Gh, 12.38. Bx 
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I. INTRODUCTION 



In the present work, we calculate the additive renormalization of the fermion mass in 
lattice QCD, using clover fermions and Symanzik improved gluons. The calculation is carried 
out up to two loops in perturbation theory and it is directly related to the determination of 
the critical value of the hopping parameter, k c . 

The clover fermion action [1] (SW) succesfully reduces lattice discretization effects and 
approaches the continuum limit faster. This justifies the extensive usage of this action in 
Monte Carlo simulations in recent years. The coefficient csw appearing in this action is a 
free parameter for the current work and our results will be given as a polynomial in csw- 

Regarding gluon fields, we employ the Symanzik improved action [2] , which also aims at 
minimizing finite lattice spacing effects. For the coefficients parameterizing the Symanzik 
action, we consider several choices of values which are frequently used in the literature. 

The lattice discretization of fermions introduces some well known difficulties; demanding 
strict locality and absence of doublers leads to breaking of chiral symmetry. In order to 
recover this symmetry in the continuum limit one must set the renormalized fermion mass 
(rrifi) equal to zero. To achieve this, the mass parameter m appearing in the Lagrangian 
must approach a critical value m c , which is nonzero due to additive renormalization. 

The mass parameter m is directly related to the hopping parameter k used in simula- 
tions. Its critical value, k c , corresponds to chiral symmetry restoration: 



where a is the lattice spacing and r is the Wilson parameter. Using Eq.(l), the non- 
renormalized fermion mass is given by: 



Thus, in order to restore chiral symmetry one must consider the limit m — > m c . This fact 
points to the necessity of an evaluation of m c . 

The perturbative value of m c is also a necessary ingredient in higher-loop calculations 
of the multiplicative renormalization of operators (see, e.g., Ref. [3]). In mass independent 
schemes, such renormalizations are typically defined and calculated at zero renormalized 
mass, and this entails setting the value of the Lagrangian mass equal to m c . 

Previous studies of the hopping parameter and its critical value have appeared in the 
literature for Wilson fermions - Wilson gluons [4] and for clover fermions - Wilson gluons [5,6]. 
The procedure and notation in our work is the same as in the above references. 

Our results for k c (and consequently for the critical fermion mass) depend on the number 
of colors (N) and on the number of fermion flavors (Nf). Besides that, there is an explicit 
dependence on the clover parameter csw which, as mentioned at the beginning, is kept as a 
free parameter. On the other hand, the dependence of the results on the choice of Symanzik 
coefficients cannot be given in closed form; instead, we present it in a list of Tables and 
Figures. 

The rest of the paper is organized as follows: In Sec. II we formulate the problem, define 
the discretized actions, and describe our calculation of the necessary Feynman diagrams. 



1 



(1) 



2 m c a + 8 r 




(2) 



2 



Sec. Ill is a presentation of our results. Finally, in Sec. IV we apply to our one- and two- 
loop results an improvement method, proposed by us [7-9]. This method resums a certain 
infinite class of subdiagrams, to all orders in perturbation theory, leading to an improved 
perturbative expansion. We end this section with a comparison of perturbative and non- 
perturbative results. Our findings are summarized in Sec. V. 



II. FORMULATION OF THE PROBLEM 

We begin with the Wilson formulation of the QCD action on the lattice, with Nf flavors 
of degenerate clover (SW) [1] fermions. In standard notation, it reads: 

S L = S G + ]T ]T(4r + m )$ f (x)il) f (x) 
f x 

~n!212 $f( x ) ( r -l^) U x,x+^f{x + fi) +ip f {x + n) (r + ^)U x+ ^ x ijf(x) 

f x,n L 

+ 7 c sw 12 12 ^Pf^a^F^x^fix), (3) 

/ x,n,v 

where : = ^ (Q^ - Q Vft ) (4) 

o 

and . Qfiu U Xj x+fj,U x -\-ij, : x+^+i'Ux-\-fj,-\-i/,x+i'Ux-\-i/,x 

U X} x+uU x+U, x+v— jJJ x+v— /i, x—fjJ^x—fj,, X 

~\~ U X; X—^JJx—jX, X — fJL — vUx — ^L — V, X — uUx — V, X 

~t~ U x ,x—vUx—v,x—v+ijJJx—v+ix,x+ijJJx+ii,,x (5) 

The clover coefficient csw is treated here as a free parameter. Particular choices of values 
for csw have been determined both perturbatively [1] and non-perturbatively [10], so as to 
minimize 0(a) effects. The Wilson parameter r is set to r = 1 henceforth; / is a flavor 
index; = j v ]. Powers of the lattice spacing a have been omitted and may be 

directly reinserted by dimensional counting. 

Regarding gluons, we use the Symanzik improved gauge field action, involving Wilson 
loops with 4 and 6 links 1 : 



2 r 

S G = — [c Re Tr { 1 - £/ p l aqU ette } + Ci 12 Re Tr { 1 - £/ re ctangle } 

9 plaquette rectangle 

+C 2 X] Re Tr {1 - t/chair} + C 3 Re Tr {1 - c7 parallclogram }] (6) 

chair parallelogram 

(g is the bare coupling constant). The lowest order expansion of this action, leading to the 
gluon propagator, is 



H x 1 plaquette, 1x2 rectangle, 1x2 chair (bent rectangle), and lxlxl parallelogram wrapped 
around an elementary 3-d cube. 
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ClJv,, 



where £ is the gauge fixing parameter (see Eq.(lO) ) and: 

G^ u (k) = kjz v + ^ {k 2 p ^nv ~ kjv p 5py\ dp P , hp = — sin — 
p a z 

= (1 - W [Co - Cx a 2 P - C 2 a 2 (kl + 
The coefficients C; are related to q by 

Co = c + 8ci + 16c 2 + 8c 3 , Ci = c 2 + c 3 , C 2 = Ci - c 2 - c 3 



(7) 



M 



(8) 



(9) 



The Symanzik coefficients must satisfy: c + 8ci + 16c 2 + 8c 3 = 1, in order to reach the 
correct classical continuum limit. Aside from this requirement, the values of q can be 
chosen arbitrarily; they are normally tuned in a way as to ensure 0(a) improvement. 

As always in perturbation theory, we must introduce an appropriate gauge-fixing term 
to the action; in terms of the gauge field Qp(x) [U XjX+ll = exp(i gQ^(x))], it reads: 



1-e 



£ Tr {A;Qp(x)A-Q u (x)}, A;Q u (x) = Q u (x -fx)- Q u (x). 



(10) 



x,p,u 



Having to compute a gauge invariant quantity, we can, for convenience, choose to work 
either in the Feynman gauge (£ = 0) or in the Landau gauge (£ = 1). Covariant gauge fixing 
produces the following action for the ghost fields uj and u7 

—^g 2 [Qp(x),[Qp{x),A+u(x)]] +•••)}, A+ou(x)=ou(x + jl)-u;(x). (11) 

Finally, the change of integration variables from links to vector fields yields a Jacobian that 
can be rewritten as the usual measure term S m in the action: 



Sm = ^Ng 2 £ £ Tr {Q,(:r)Q M (:r)} + 



(12) 



In and 5 m we have written out only terms relevant to our computation. The full action 
is: S = Sl + S g f + Sgh + S m . 

The bare fermion mass m# must be set to zero for chiral invariance in the classical 
continuum limit. Terms proportional to r in the action, as well as the clover terms, break 
chiral invariance. They vanish in the classical continuum limit; at the quantum level, they 
induce nonvanishing, flavor-independent fermion mass corrections. Numerical simulation 
algorithms usually employ the hopping parameter, 
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K = 



2 m a + 8 r 
4 



(13) 



as an adjustable input. Its critical value, at which chiral symmetry is restored, is thus l/8r 
classically, but gets shifted by quantum effects. 

The renormalized mass can be calculated in textbook fashion from the fermion self- 
energy. Denoting by S L (jo, m ,g) the truncated, one particle irreducible fermion two-point 
function, we have for the fermion propagator: 

S(p) = [i$ +m(p) - i: L (p,m ,g) (14) 

1 2r 
where : $ = — J^7 M sin(ap M ), m(p) = m Q H V] sin 2 (ap M /2) 

a ij, a n 

To restore the explicit breaking of chiral invariance, we require that the renormalized 
mass vanish: 

5_1 (°)| mo ^ mc = =► m c = J: L (0,m c ,g) (15) 

The above is a recursive equation for m c , which can be solved order by order in perturbation 
theory. 

We denote by dm the additive mass renormalization of m : m B = m Q — dm. In or- 
der to obtain a zero renormalized mass, we must require m B — > 0, and thus m Q — * dm. 
Consequently, 

m c = dm = dm(i_i oop ) + dm^-ioo P ) (16) 

At tree level, m c = 0. 

Two diagrams contribute to rfm(i_i oop ), shown in Fig. 1. In these diagrams, the fermion 
mass must be set to its tree level value, m Q — > 0. 




Fig. 1. One- loop diag rams contributing to d?n^_ i 00 p). 
Wavy (solid) lines represent gluons (fermions). 

The quantity dm^-ioop) receives contributions from a total of 26 diagrams, shown in Fig. 
2. Genuine two-loop diagrams must again be evaluated at m — > 0; in addition, one must 
include to this order the one-loop diagram containing an 0(g 2 ) mass counterterm (diagram 
23). 

Certain sets of diagrams, corresponding to one-loop renormalization of propagators, must 
be evaluated together in order to obtain an infrared convergent result: These are diagrams 
7+8+9+10+11, 12+13, 14+15+16+17+18, 19+20, 21+22+23. 
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Fig. 2. Two-loop diag rams contributing to dm^— loop) • 

Wavy (solid, dotted) lines represent gluons 
(fermions, ghosts). Crosses denote vertices stemming from the measure part of the action; a solid 
circle is a fermion mass counterterm. 

III. COMPUTATION AND RESULTS 

Given that the dependence of m c on the Symanzik coefficients q cannot be expressed 
in closed form, we chose certain sets of values for q , presented in Table I, which are in 
common use [11-16]: Plaquette, Symanzik (tree level improved), Tadpole Improved Liischer- 
Weisz (TILW), Iwasaki and DBW2. Actually, since the gluon propagator contains only the 
combinations C\ and C 2 (Eq.(9)), all results for m c can be recast in terms of C 1: C 2 and 
one additional parameter, say, c 2 ; in this case the dependence on c 2 (at fixed Ci, C 2 ) is 
polynomial of second degree. 

The contribution drrii of the I th one-loop diagram to dm, can be expressed as: 
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(^ 2 - h r , j „ w 

.(0 



^ = ^V^^-E4w^ (17) 



i=0 



where e\ are numerical one-loop integrals whose values depend on Ci, C2. The dependence 
on csw is seen to be polynomial of degree 2 (i — 0, 1, 2). 

The contribution to (im from two-loop diagrams that do not contain closed fermion loops, 
can be written in the form 

dm = g 4 • £ 4w N* 4 ef ^ (18) 

where the index / runs over all contributing diagrams, j = 0, 2 and = 0, 1, 2 (since up to 
two vertices from the gluon action may be present in a Feynman diagram). The dependence 
on csw is now polynomial of degree 4 (i = 0, 4). The coefficients e^' -7 ''*^ (as well as 
ef^ of Eq.(19) below) are two-loop numerical integrals; once again, they depend on Ci, C 2 . 
Finally, the contribution to dm from two-loop diagrams containing a closed fermion loop, 
can be expressed as 

^ = ^^A^ 4 -£c* sw ef (19) 
JV i=0 

where the index I runs over diagrams 12-13, 19-20. Summing up the contributions of all 
diagrams, dm assumes the form 

dm = J2 dm, = , 2 ■ £ c* sw e W + , 4 • £ 4w ^ 4 e ( - fc) 

i ' i ' i,j,k 
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+^V^ Ar ^ 4 -^ ciswg(i) (2o) 

i 

In the above ; e^ ,k \ are the sums over all contributing diagrams of the quantities: 
ef\ ef' i,fe) , ef\ respectively (cf. Eqs.(17,18,19) ). 

The coefficients lead to the total contribution of one-loop diagrams. Their values are 
listed in Table II, for the ten sets of q values shown in Table I. Similarly, results for the 
coefficients e^'^ and corresponding to the total contribution of two-loop diagrams, are 
presented in Tables III- VII. 

In order to enable cross-checks and comparisons, numerical per-diagram values of the 
constants ef\ ef^'^ and ejp are presented in Tables VIII-XII, for the case of the Iwasaki 
action. For economy of space, several vanishing contributions to these constants have simply 
been omitted. A similar breakdown for other actions can be obtained from the authors upon 
request. 

The total contribution of one-loop diagrams, for N = 3 can be written as a function of 
the clover parameter csw- In the case of the Plaquette, Iwasaki, and DBW2 actions, we 
find, respectively: 

rfm (i-Zp) e = 9 2 ( - 0.434285489(1) + 0.1159547570(3) c sw + 0.0482553833(1) c 2 sw ) (21) 
dmfn^p) = g 2 ( - 0.2201449497(1) + 0.0761203698(3) c sw + 0.0262264231(1) c 2 sw ) (22) 
rfm fi-^o P ) = 9 2 ( - 0.0972070995(5) + 0.0421775310(1) c sw + 0.01141359801(1) cg w ) (23) 
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A similar process can be followed for two-loop diagrams. In this case, we set N = 3, 
C2 = and we use three different values for the flavor number: Nf = 0, 2, 3. Thus, for the 
Plaquette, Iwasaki and DBW2 actions, the total contribution is, respectively: 



N f = 0: 



N f 



2 : 



Nt = 3 : 



7 Plaquette 
am (2-loop) 



* 4 ( 
» 4 ( 



0.1255626(2) + 0.0203001(2) c sw + 0.00108420(7) c 2 sw 
0.00116538(2) 4 W - 0.0000996725(1) c| w ) (24) 
0.1192361(2) + 0.0173870(2) c sw + 0.00836498(8) c 2 sw 
0.00485727(3) c| w - 0.001 1561947(4) c| w ) (25) 
0.1160729(2) + 0.0159305(2) c sw + 0.0120054(1) c 2 sw 
0.00670321(3) cg w - 0.0016844558(6) c| w ) (26) 



iV / = 0: 



JV/ = 2 : 



iV/ = 3 : 



^ Iwasaki 
am (2-loop) 



Am Iwasaki 
« m (2-loop) 



j Iwasaki 
« m (2-loop) 



* 4 ( 



0.0099523(2) - 0.0024304(5) c sw - 0.00232855(4) c 2 sw 
0.00032100(2) 4 W - 0.0000419365(1) c 4 sw ) (27) 
0.0076299(2) - 0.0040731(5) c sw + 0.00102758(6) c 2 sw 
0.00242924(3) c 3 sw - 0.000457690(2) c| w ) (28) 
0.0064687(2) - 0.0048944(5) c sw + 0.00270565(7) c 2 sw 
0.00348335(3) 4 W - 0.000665567(2) c| w ) (29) 



.2 

-sw 



N f = : dmg^op) = g 4 ( + 0.005099(2) - 0.0053903(7) c sw - 0.0011157(1) 

- 0.00004482(2) c 3 sw - 0.0000111470(2) c 4 sw ) (30) 
rf"i(2-to?p) = ^ 4 ( + 0.005944(2) - 0.0061840(7) c sw + 0.0002046(2) c| w 

- 0.0010177(3) 4 W - 0.000125065(3) c^ w ) (31) 
^f 2 -too P ) = 9 4 ( + 0.006366(2) - 0.0065809(7) c sw + 0.0008648(2) c^ w 

- 0.0015042(4) 4 W - 0.000182023(5) c 4 sw ) (32) 



N f = 2 



N f = 3 : 



In Figs. 3, 4, and 5 we present the values of dm(2-\ oop ) for Nf = 0, 2, 3, respectively; 
the results are shown for all choices of Symanzik actions which we have considered, as a 
function of Csw (N = 3, c 2 = 0). In all cases, the dependence on csw is rather mild. One 
observes that dm^-ioop) is significantly smaller for all improved actions, as compared to the 
plaquette action; in particular, in the case of DBW2, dm^-ioop) is closest to zero and it 
vanishes exactly around csw = 1- 

Another feature of these results is that they change only slightly with Nf , especially in 
the range csw < 1-5 . This is due to the small contributions of diagrams with closed fermion 
loops (diagrams 12, 13, 19, 20). By the same token, in the case of nondegenerate flavors, 
dm(2-i 00 p) is expected to depend only weakly on the mass of the virtual fermion. 
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-sw 
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Fig. 3. Total contribution of two-loop diagrams, for N = 3, Nf = and C2 = 0. Legends appear 
in the same top-to-bottom order as the corresponding lines. 



0.02 




"SW 



Fig. 4. Total contribution of two-loop diagrams, for N = 3, Nf = 2 and C2 = 0. Legends appear 
in the same top-to-bottom order as the corresponding lines. 
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0.02 




-0.12 1 1 ' 1 1 

0.0 0.5 1.0 1.5 2.0 



c sw 

Fig. 5. Total contribution of two-loop diagrams, for N = 3, Nf = 3 and C2 = 0. Legends appear 
in the same top-to-bottom order as the corresponding lines. 

IV. IMPROVED PERTURBATION THEORY 

We now apply our method of improving perturbation theory [7-9], based on resummation 
of an infinite subset of tadpole diagrams, termed "cactus" diagrams. In Ref. [9] we show 
how this procedure can be applied to any action of the type we are considering here, and 
it provides a simple, gauge invariant way of dressing, to all orders, perturbative results at 
any given order (such as the one- and two-loop results of the present calculation). Some 
alternative ways of improving perturbation theory have been proposed in Refs. [17,18]. In 
a nutshell, our procedure involves replacing the original values of the Symanzik and clover 
coefficients by improved values, which are explicitly computed in [9]. Applying at first this 
method to one-loop diagrams, the improved ("dressed") value dm dr of the critical mass 
(N = 3, C2 = 0) can be written as: 

2 

^ m (l-loop) = X! E dr C SW (33) 
i=0 

In comparing with eW of Eq. (20), the quantity s d % l is the result of one-loop Feynman dia- 
grams with dressed values for the Symanzik parameters, and it has already been multiplied 
by g 2 (N 2 — 1) /N. The dependence of 

£ dr on 9 is quite complicated now, and cannot be 
given in closed form; instead e^. must be computed numerically for particular choices of g. 
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Listed in Table XIII are the results for e$l along with the value of (3 — 2N/g 2 corresponding 
to each one of the 16 actions used in this calculation. 

An attractive feature of this improvement procedure is that it can be applied also to 
higher loop perturbative results, with due care to avoid double counting of the cactus dia- 
grams which were already included at one loop. Ideally, of course, one loop improvement 
should already be adequate enough, so as to obviate the need to consider higher loops; 
indeed, we find this to be the case and, consequently, we limit our discussion of two-loop 
improvement to only the plaquette action (/3 = 5.29, N = 3, Nf = 2), the Iwasaki action 
((3 = 1.95, iV = 3, N f = 2) and the DBW2 action ((3 = 0.87 and (3 = 1.04, N = 3, N f = 2). 
Using these values, the contribution to cfrft^-ioop) * s a polynomial in csw : 

<H r -ioo P ), plaquette = "0.77398(8) + 0.16330(4) c sw + 0.06224534(1) c 2 sw 

-0.0044006(9) 4 W - 0.00073780(6) c^ w (34) 

^(2-ioop), iwasaki = -0.0813302(9) + 0.043030(3) c sw + 0.0308196(2) c 2 sw 

-0.00767090(8) c| w - 0.001160923(1) 4 W (35) 

^( 2 r -ioop),DBW2 {/3 =o.87) = -0.044906(1) + 0.029449(4) c sw + 0.0239522(2) c 2 sw 

-0.0082231(1) 4 W - 0.001218955(4) c 4 sw (36) 

^(2-ioop),DBW2(/3=i.o4) = -0.031260(1) + 0.021793(2) c sw + 0.0188027(2) c 2 sw 

-0.00705284(9) c| w - 0.001055657(1) c| w (37) 

The comparison between the total dressed contribution dm Ar = dm9l_ loop ^ + ^^(2-ioop) 
and the unimproved contribution, dm, for the plaquette action is exhibited in Fig. 6, as a 
function of csw- Similarly, dm dr for the Iwasaki and the DBW2 actions is shown in Fig. 7 
and Fig. 8, respectively. 
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c sw 

Fig. 8. Improved and unimproved values of dm up to two loops, as a function of csw; f° r the 
DBW2 action (N = 3, JV> = 2). We set j3 = 0.87 (solid lines) and (3 = 1.04 (dotted lines). 

Finally, in Table XIV, we present a comparison of dressed and undressed results, for some 
commonly used values of /3, Nf, csw ; and we also compare with available non perturbative 
estimates for k c [10,19-22]. We observe that improved perturbation theory, applied to one- 
loop results, already leads to a much better agreement with the non perturbative estimates. 

V. DISCUSSION 

To recapitulate, in this paper we have calculated the critical mass m c , and the associated 
critical hopping parameter k c , up to two loops in perturbation theory, using the clover 
action for fermions and the Symanzik improved gluon action with 4- and 6-link loops. The 
perturbative value of m c is a necessary ingredient in the higher-loop renormalization of 
operators, in mass independent schemes: Such renormalizations are typically defined and 
calculated at vanishing renormalized mass, which amounts to setting the Lagrangian mass 
equal to m c . 

In our calculations, we have chosen for the Symanzik coefficients q a wide range of 
values, which are most commonly used in numerical simulations. The dependence of our 
results on the number of colors N and the number of fermion flavors Nf is shown explicitly. 
The dependence on the clover parameter csw is hi the form of a fourth degree polynomial 
whose coefficients we compute explicitly; it is expected, of course, that the most relevant 
values for csw are those optimized for 0(a) improvement, either at tree level (csw = 1), or 
at one loop [1], or non-perturbatively [10]. 

Since m c is gauge invariant, we chose to calculate it in the Feynman gauge. The propaga- 
tor appearing in Feynman diagrams is the inverse of a nondiagonal matrix; while this inverse 
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can be written down explicitly, it is more convenient, and more efficient in terms of CPU 
time, to perform the inversion numerically. Integrations over loop momenta were performed 
as momentum sums on lattices of finite size L, where typically L < 40; extrapolation to 
L — > oo introduces a systematic error, which we estimate quite accurately. 

Our results for m c are significantly closer to zero in the case of Symanzik improved 
actions, as compared to the plaquette action. In particular, the DBW2 action stands out 
among the rest, in that m c vanishes exactly for a value of csw around 1. Thus, improved 
actions seem to bring us quite near the point of chiral symmetry restoration. The dependence 
of m c on the number of flavors is seen to be very mild. This fact would also suggest that, in 
the case of nondegenerate flavors, m c should depend only weakly on the mass of the virtual 
fermion. 

Finally, we have made some comparisons among perturbative and non-perturbative re- 
sults for k c . While these are expected to differ for a power divergent additive renormaliza- 
tion, such as the quantity under study, we nevertheless find a reasonable agreement. This 
agreement is further enhanced upon using an improved perturbative scheme, which entails 
resumming, to all orders in the coupling constant, a dominant subclass of tadpole diagrams. 
The method, originally proposed for the Plaquette action (see Ref. [7]), was extended in 
Ref. [9] to encompass all possible gluon actions made of closed Wilson loops, and can be 
applied at any given order in perturbation theory. As would be desirable, one-loop improve- 
ment is seen to be already adequate to give a reasonable agreement among perturbative and 
non-perturbative values. Indeed, our results for nfL\ oop are significally closer to the non- 
perturbative evaluations, as shown in Table XIV; in fact, the two- loop dressing procedure 
introduces no further improvement to the comparison. 

Acknowledgments: This work is supported in part by the Research Promotion Founda- 
tion of Cyprus (Proposal Nr: ENTAS/0504/11, ENIEX/0505/45). 
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TABLES 





TABLE I. Input parameters cq, 


Cl, 


c 3 (C2 = 0) 






Action 




c 




Ci 






Plaquette 




1.0 


n 
u 




n 
u 




Symanzik 




1.6666667 


-0. 


.083333 







TILW, /3c = 


8.60 


2.3168064 


-0. 


151791 


-0. 


0128098 


TILW, /3c = 


8.45 


2.3460240 


-0. 


.154846 


-0. 


0134070 


TILW, (3c = 


8.30 


2.3869776 


-0. 


.159128 


-0. 


0142442 


TILW, (3c = 


8.20 


2.4127840 


-0. 


.161827 


-0. 


0147710 


TILW, (3c = 


8.10 


2.4465400 


-0. 


.165353 


-0. 


.0154645 


TILW, (3c = 


8.00 


2.4891712 


-0. 


.169805 


-0. 


0163414 


Iwasaki 




3.648 


-0. 


.331 







DBW2 




12.2688 


-1. 


.4086 








TABLE II. Total contribution of one-loop diagrams 



Action 



Plaquette 
Symanzik 
TILW (8.60) 
TILW (8.45) 
TILW (8.30) 
TILW (8.20) 
TILW (8.10) 
TILW (8.00) 
Iwasaki 
DBW2 



,(o) 



-0.1628570582(5) 

-0.12805490528(8) 

-0.10821568768(4) 

-0.10749185625(3) 

-0.1064962872(3) 

-0.1058799831(2) 

-0.1050866191(1) 

-0.10410447893(3) 

-0.08255435613(4) 

-0.0364526623(2) 



0.0434830339(1) 

0.0378314931(2) 

0.03408560232(6) 

0.0339409375(1) 

0.0337409869(2) 

0.0336166372(1) 

0.03345591621(5) 

0.03325593631(8) 

0.0285451387(1) 

0.01581657412(5) 



r(2) 



0.01809576875(4) 

0.01476335801(5) 

0.01265991972(4) 

0.01258108895(1) 

0.012472434543(4) 

0.0124050416(1) 

0.012318127134(5) 

0.012210297749(7) 

0.00983490867(5) 

0.004280099253(2) 
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TABLE III. Total contribution of two-loop diagrams of order 0(N 2 , c[>) 



Action 

Plaquette 
Symanzik 
TILW (8.60) 
TILW (8.45) 
TILW (8.30) 
TILW (8.20) 
TILW (8.10) 
TILW (8.00) 
Iwasaki 
DBW2 



e (0,2,0) 

-0.01753602(2 
-0.00810366(1 
-0.00437013(7 
-0.00425575(7 
-0.00410086(7 
-0.00400636(6 
-0.00388630(6 
-0.00374009(6 
-0.00112957(2 
0.0008481(2) 



e (l,2,0) 

0.00259963(2 
0.00095046(2 
0.00019403(5 
0.00016978(6 
0.00013682(7 
0.00011666(8 
0.00009097(9 
0.00005958(9 
-0.00052964(6 
-0.00085301(8 



e (2,2,0) 

-0.000155894(8) 

-0.000404510(9) 

-0.00045894(1) 

-0.00045962(1) 

-0.00046040(1) 

-0.00046080(1) 

-0.00046123(1) 

-0.000461601(9) 

-0.000436966(5) 

-0.00018540(1) 



e (3,2,0) 

-0.000163242(2) 
-0.000107348(2) 
-0.000078117(3) 
-0.000077102(3) 
-0.000075713(3) 
-0.000074857(3) 
-0.000073760(3) 
-0.000072410(3) 
-0.000045009(3) 
-0.000006164(3) 



e (4,2,0) 

-0.00001721759(2 
-0.00001275904(1 
-0.00001020820(1 
-0.00001011451(1 
-0.00000998564(1 
-0.00000990584(1 
-0.00000980314(1 
-0.00000967600(1 
-0.00000682353(1 
-0.00000173502(3 



TABLE IV. Total contribution of two-loop diagrams of order O(N , d]) 



Action 

Plaquette 
Symanzik 
TILW (8.60) 
TILW (8.45) 
TILW (8.30) 
TILW (8.20) 
TILW (8.10) 
TILW (8.00) 
Iwasaki 
DBW2 



e (o,o,o) 

0.01656633(2 
0.00605656(1 
0.00202637(3 
0.00190729(3 
0.00174666(3 
0.00164901(3 
0.00152532(3 
0.00137535(4 
-0.00103022(1 
-0.0018961(2) 



e (i,o,o) 

-0.00055904(1) 
0.000935801(6) 
0.00157890(3) 
0.00159800(3) 
0.00162375(2) 
0.00163939(2) 
0.00165917(2) 
0.00168310(3) 
0.00203254(1) 
0.0016130(3) 



e (2,0,0) 

0.002622771(7 
0.002120980(9 
0.001790242(9 
0.001777415(9 
0.001759689(9 
0.001748661(9 
0.001734421(9 
0.00171671(1) 
0.001313076(3 
0.000413397(9 



e (3,0,0) 

0.000158125(2 
0.000104973(2 
0.000076167(2 
0.000075164(3 
0.000073791(3 
0.000072944(3 
0.000071859(3 
0.000070522(3 
0.000043949(3 
0.000005057(3 



e (4,o,o) 

0.00004282674(2 
0.00002971553(1 
0.00002260669(1 
0.00002235603(1 
0.00002201243(1 
0.00002180041(1 
0.00002152826(1 
0.00002119259(1 
0.00001423324(1 
0.00000307480(3 



TABLE V. Total contribution of two-loop diagrams containing closed fermion loops 



Action 



5(0) 



5(1) 



;(2) 



2(3) 



5(4) 



Plaquette 

Symanzik 

TILW (8.60) 

TILW (8.45) 

TILW (8.30) 

TILW (8.20) 

TILW (8.10) 

TILW (8.00) 

Iwasaki 

DBW2 



.00118621(2 
00081496(1 
00063643(1 
00063033(1 
00062198(1 
00061684(1 
00061025(1 
00060214(1 
00043546(1 
.00015833(3 



-0.000546197(8) 

-0.000448276(6) 

-0.000389464(5) 

-0.000387269(5) 

-0.000384243(5) 

-0.000382366(5) 

-0.000379946(5) 

-0.000376945(5) 

-0.00030800(1) 

-0.00014883(4) 



0.001365146(9 
0.001041379(8 
0.000857737(3 
0.000851127(3 
0.000842047(3 
0.000836433(3 
0.000829214(4 
0.000820289(4 
0.000629274(8 
0.00024756(2) 



-0.000692228(3) 
-0.000574521(3) 
-0.000500011(5) 
-0.000497194(5) 
-0.000493307(5) 
-0.000490894(5) 
-0.000487781(4) 
-0.000483915(4) 
-0.000395294(3) 
-0.00018242(5) 



-0.00019809791(7) 

-0.0001453370(2) 

-0.0001148491(1) 

-0.0001137544(1) 

-0.0001122515(1) 

-0.0001113227(1) 

-0.0001101288(1) 

-0.0001086536(1) 

-0.0000779538(3) 

-0.0000213595(6) 
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TABLE VI. Total contribution of two-loop diagrams containing the parameter C2 (part 1) 



Action 

Plaquette 
Symanzik 
TILW (8.60) 
TILW (8.45) 
TILW (8.30) 
TILW (8.20) 
TILW (8.10) 
TILW (8.00) 
Iwasaki 
DBW2 



e (o,o,i) 

0.077167(3) 

0.034929(2) 

0.020247(1) 

0.019816(1) 

0.019235(1) 

0.018881(1) 

0.018433(1) 

0.017888(1) 

0.0087615(7) 

0.0007907(2) 



e (i,o,i) 

-0.019808(3) 
-0.010895(2) 
-0.007117(2) 
-0.006998(2) 
-0.006835(2) 
-0.006736(2) 
-0.006609(2) 
-0.006454(2) 
-0.003656(1) 
-0.0004889(3) 



e (2,0,l) 

-0.0085415(2) 

-0.0041454(2) 

-0.0024559(1) 

-0.0024050(1) 

-0.0023362(1) 

-0.0022942(1) 

-0.0022410(1) 

-0.0021762(1) 

-0.00107856(8) 

-0.00008343(2) 



e (0,2,l) 

-0.047102(4) 
-0.017940(2) 
-0.008702(1) 
-0.008448(1) 
-0.0081078(6) 
-0.0079023(7) 
-0.0076431(9) 
-0.0073300(6) 
-0.0027484(4) 
0.0001308(2) 



e (1.2,l) 

0.010439(3) 
0.004491(2) 
0.002251(1) 
0.002185(1) 
0.0020973(9) 
0.002044(1) 
0.0019761(8) 
0.0018940(6) 
0.0006646(5) 
-0.0001587(3) 



TABLE VII. Total contribution of two-loop diagrams containing the parameter C2 (part 2) 



Action 

Plaquette 
Symanzik 
TILW (8.60) 
TILW (8.45) 
TILW (8.30) 
TILW (8.20) 
TILW (8.10) 
TILW (8.00) 
Iwasaki 
DBW2 



e (2,2,l) 

0.0039245(3) 
0.0014622(1) 
0.0006472(1) 
0.0006251(1) 
0.0005954(1) 
0.0005775(1) 
0.0005550(1) 
0.0005279(1) 
0.00015719(6) 
-0.00002436(1) 



e (3,2,l) 

-0.0000842143(1) 
-0.0000454986(1) 
-0.00002872341(6 
-0.00002818123(6 
-0.00002744385(5 
-0.00002699223(5 
-0.00002641646(5 
-0.00002571231(5 
-0.00001249281(2 
-0.00000050404(9 



e (0,2,2) 

-0.09448252(9) 

-0.03417549(2) 

-0.017374635(6) 

-0.016917713(6) 

-0.016304614(5) 

-0.015933835(5) 

-0.015466270(5) 

-0.014902324(4) 

-0.00596123(2) 

-0.00028731(2) 



e (l,2,2) 

0.02755993(3) 

0.01248953(1) 

0.007205477(3) 

0.007049188(2) 

0.006838088(3) 

0.006709626(4) 

0.006546741(4) 

0.006348924(5) 

0.00295502(1) 

0.00020317(4) 



e (2,2,2) 

0.010521016(1) 
0.0041047891(2 
0.0021218443(2 
0.0020666192(2 
0.0019924047(3 
0.0019474604(2 
0.0018907121(3 
0.0018221643(3 
0.0007286816(4 
0.0000278810(8 



TABLE VIII. Contribution of one-loop diagrams, for the Iwasaki action 

~. (ij (i) 

* £i £2 

-0.05602636832(2) -0.02652798781(3) 

1 0.0285451387(1) 

2 0.00983490867(5) 



TABLE IX. Contribution of diagrams 3, 4, 6, for the Iwasaki action 



k 



e 4 



-0.0003923686(9) 

2 0.0002615791(6) 

10 

12 

2 

2 2 



-0.000743134(3) 
0.000495422(2) 
0.001900337(2) 
0.0017774410(9) 
-0.0010339720(2) 
-0.001041123(1) 



-0.0000714882(8) 
0.0000357441(4) 




0.0002799238(4) 
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TABLE X. Contribution of diagrams 7-11, 14-18, 24, 26, for the Iwasaki action 



i 


3 


K 


e 7-ll 


e 14-18 


e 24 


e 26 











0.00042802(1) 


-0.000195263(2) 














1 


0.0057103(7) 


0.0030512(2) 











2 





-0.00111995(2) 


-0.00029748(1) 





-0.000298742(2) 





2 


1 


-0.0022472(3) 


-0.0008718(2) 





0.0003705893(7) 





2 


2 


-0.00371263(2) 


-0.00224859(1) 








1 











0.00064534(1) 








1 





1 





-0.003656(1) 








1 


2 








0.00011079(6) 


-0.000144897(2) 


0.000429899(1) 


1 


2 


1 





0.0006450(5) 


0.000248682(4) 


-0.00022905(1) 


1 


2 


2 





0.00295502(1) 








2 











-0.000000974(1) 








2 





1 





-0.00107856(8) 








2 


2 








0.000141960(3) 


0.000042314(2) 


0.0003303085(7) 


2 


2 


1 





0.00039546(6) 


0.00002909398(7) 


-0.000267364(2) 


2 


2 


2 





0.0007286816(4) 








3 


2 














-0.000019835(1) 


3 


2 


1 











-0.00001249281(2) 



TABLE XI. Contribution of diagrams 12, 13, 19, 20, for the Iwasaki action 



i 


g (<) 


e 19-20 





0.000261920(6) 


0.000173538(9) 


1 


-0.0000308339(1) 


-0.00027717(1) 


2 


0.000370942(2) 


0.000258332(8) 


3 





-0.000395294(3) 


4 





-0.0000779538(3) 



TABLE XII. Contribution of diagrams 21-23, 25, 27, 28, for the Iwasaki action 



i 


3 


k 


e 21-23 


_(»j>fc) 

e 25 


e 27 


e 28 











0.000373419(3) 


-0.000158621(4) 


-0.000094848(3) 


-0.0001759336(5) 





2 





-0.000373419(3) 


0.000079311(2) 





0.0000879668(3) 


1 








-0.000887295(1) 


0.0001396819(4) 


0.000045158(4) 


0.0001893113(5) 


1 


2 





0.000887295(1) 


0.000085189(2) 





-0.000120480(1) 


2 








0.000194437(1) 


-0.0000319392(3) 


0.000168506(2) 


-0.0000509266(2) 


2 


2 





-0.000194437(1) 


-0.000005787(2) 





0.0000098758(1) 


3 








0.000059183(3) 





-0.000015234(1) 





3 


2 





-0.000059183(3) 


0.0000172022(5) 





0.0000168072(6) 


4 








0.00000682353(1) 





0.000007409712(6) 





4 


2 





-0.00000682353(1) 
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TABLE XIII. Results for dmf[_ lo } (Eq. (33)), with N = 3 



Action 


& 


e (0) 

^dr- 
ill 


e {1) 

C dr- 
ill 


e (2) 

Ltl 


Plaquette 


6.00 


-0.579221119(2) 


0.1159547570(3) 


0.03618067788(9) 


Symanzik 


5.00 


-0.4869797578(8) 


0.1121369999(4) 


0.03538605357(4) 


Symanzik 


5.07 


-0.478756110(2) 


0.11072412996(5) 


0.03507238306(5) 


Symanzik 


6.00 


-0.3915226522(2) 


0.0947962001(5) 


0.03124138429(9) 


TILW (8.60) 


3.7120 


-0.5358770348(7) 


0.1265917638(3) 


0.03813963851(4) 


TILW (8.45) 


3.6018 


-0.5497415338(3) 


0.1291104644(3) 


0.0386337113(1) 


TILW (8.30) 


3.4772 


-0.5651407386(9) 


0.1319263769(1) 


0.0391695069(1) 


TILW (8.20) 


3.3985 


-0.5756111531(9) 


0.1337937558(7) 


0.03951713046(7) 


TILW (8.10) 


3.3107 


-0.5870122772(4) 


0.1358437825(6) 


0.0398899143(3) 


TILW (8.00) 


3.2139 


-0.599415804(1) 


0.138085996(2) 


0.0402877133(4) 


Iwasaki 


1.95 


-0.757856451(1) 


0.1671007819(8) 


0.044746728234(1) 


Iwasaki 


2.20 


-0.6555102085(5) 


0.1537748193(6) 


0.04293183656(3) 


Iwasaki 


2.60 


-0.541348980(1) 


0.1359882440(3) 


0.03967626495(6) 


DBW2 


0.6508 


-0.7749943512(7) 


0.1847244889(1) 


0.04731717866(3) 


DBW2 


0.8700 


-0.574781578(1) 


0.1575688409(9) 


0.04281261980(1) 


DBW2 


1.0400 


-0.4822863343(9) 


0.1412499230(5) 


0.039186543574(5) 



TABLE XIV. One- and two-loop results, and non-perturbative estimates for k c 



Action 


N f 




csw 


«i 


—loop 


K 5 


!— loop 


dr 

K l-loop 


K 2-loop 


„non- pert 


[Ref.] 


Plaquette 





6.00 


1 


.479 


0. 


1301 


0. 


1335 


0. 


1362 


0, 


1362 


0.1392 


[19] 


Plaquette 





6.00 


1 


.769 


0. 


1275 


0. 


1306 


0. 


1337 


0. 


1332 


0.1352 


[10] 


Plaquette 


2 


5.29 


1 


.9192 


0. 


1262 


0. 


1307 


0. 


1353 


0. 


1341 


0.1373 
0.1363 


[20] 
[21] 


Iwasaki 


2 


1.95 


1 


.53 


0. 


1292 


0. 


1368 


0. 


1388 


0. 


1379 


0.1421 


[22] 


TILW (8.60) 





3.7120 


1 


.0 


0. 


1339 


0. 


1370 


0. 


1378 


0. 


1384 






TILW (8.00) 





3.2139 


1 


.0 


0. 


1348 


0. 


1387 


0. 


1397 


0. 


1406 






DBW2 


2 


0.87 





.0 


0. 


1502 


0. 


1384 


0. 


1460 


0. 


1479 






DBW2 


2 


0.87 


1 


.0 


0. 


1352 


0. 


1372 


0. 


1379 


0. 


1379 






DBW2 


2 


1.04 





.0 


0. 


1454 


0. 


1375 


0. 


1421 


0. 


1434 






DBW2 


2 


1.04 


1 


.0 


0. 


1334 


0. 


1348 


0. 


1352 


0. 


1352 
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